Secondary atmospheric tau neutrino production 
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We evaluate the flux of tau neutrinos produced from the decay of pair produced taus from incident 
muons using a cascade equation analysis. To solve the cascade equations, our numerical result for the 
tau production Z moment is given. Our results for the flux of tau neutrinos produced from incident 
muons are compared to the flux of tau neutrinos produced via oscillations and the direct prompt 
atmospheric tau neutrino flux. Results are given for both downward and upward going neutrinos 
fluxes and higher zenith angles are discussed. We conclude that the direct prompt atmospheric tau 
neutrino flux dominates these other atmospheric sources of tau neutrinos for neutrino energies larger 
than a few TeV for upward fluxes, and over a wider range of energy for downward fluxes. 
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I. INTRODUCTION 

Atmospheric muons and muon neutrinos dominate the 
atmospheric lepton flux for a lepton energy of 100 GeV 
[l|. For E = 100 GeV, the dominant sources of atmo- 
spheric leptons are pion and kaon production by cos- 
mic ray interactions in the atmosphere (the "conventional 
flux") [2] at a height of 15-20 km 0, followed by their 
decay. Cosmic ray production of charm is suppressed at 
this energy, where the "prompt flux" of muon neutrinos 
is about three orders of magnitude smaller than the con- 
ventional flux of muon neutrinos [J, |5|. Charm decays, 
via the Dg meson, can also produce tau neutrinos. The 
prompt tau neutrino flux is about one order of magnitude 
smaller than the prompt muon neutrino flux [6J, |7| . 

Neutrino oscillations of atmospheric neutrinos can pro- 
duce appreciable atmospheric tau neutrino fluxes in some 
energy ranges. Indeed, the oscillations of upward muon 
neutrinos to tau neutrinos as they propagate along the 
diameter of the Earth resulted in a convincing deficit 
in muon neutrinos in the sub-GeV and multi-GeV event 
samples in the SuperKamiokande detector Q. Down- 
ward atmospheric neutrinos have a relatively short dis- 
tance to propagate, so the oscillation probability is re- 
duced. At higher energies, neutrino oscillations are also 
suppressed. 

Neutrino oscillations over astrophysical distances will 
provide for astrophysical tau neutrino fluxes. The neutri- 
nos may come from pion and kaon production and decay 
in sources of cosmic rays, or they may come from more 
exotic sources, such as dark matter annihilation to neu- 
trinos or particles which decay to neutrinos. 

In this Brief Report, we explore another source of tau 
neutrinos, namely tau pair production by atmospheric 
muons in transit through the Earth. These tau pairs 
may lose energy as they traverse a portion of the Earth, 
but for low enough energy, they do not since the tau life- 
time is so short [9| . We provide here an evaluation of the 
tau neutrino flux from these tau pairs which we compare 
with the oscillated atmospheric and prompt atmospheric 
tau neutrino fluxes. We find that the fiux of tau neutri- 
nos from muon production of tau pairs is higher than the 
oscillated tau neutrino flux in the downward direction 



for energies above ~ 1000 GeV. In the upward direction, 
the oscillated flux is higher than that from muon produc- 
tion out to ^ 10^ GeV. In both cases, however, the tau 
neutrino flux from muon production of tau pairs is lower 
than the prompt atmospheric tau neutrino flux. 

In the next section we describe how the flux of tau 
neutrinos from tau pairs is evaluated. In Sec. Ill, we 
discuss the oscillated and prompt tau neutrino fluxes. 
Our results for the downward and upward tau neutrino 
fluxes are shown in Sec. IV. 



II. TAU NEUTRINOS FROM ATMOSPHERIC 
MUON TAU PAIR PRODUCTION 



Muon production of electron-positron pairs as muons 
pass through materials is one of the dominant energy loss 
mechanisms for muons. While kinematically suppressed, 
tau pairs are also produced. Up to Er = 10^ GeV and 
higher, tau energy loss does not compete with tau decay 
Q, with each decay resulting in a tau neutrino. While 
the tau pair production cross section is lower than the 
electron-positron pair production cross section by a fac- 
tor of'-^lO"^ — 10~* [ly at 1 TeV muon energies, there 
are orders of magnitude more muons than tau neutrinos 
coming from direct production or from oscillations in the 
downward direction. 

To evaluate the tau neutrino flux from muon produc- 
tion of tau pairs, we use the Z-moment method often ap- 
plied to atmospheric neutrino calculation [^. The start- 
ing point is the coupled set of cascade equations in terms 
of differential lepton fluxes (j)i which depend on the lepton 
energy and the column depth X, 
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Here, we use the interaction length which depends 
on the tau pair production cross section \pair{E) = 
A/{NA(Jpair{E)) [10] and the decay length Xdec{E) = 
pE/{Tr m-r) where T^ = I/tt, the inverse lifetime of the 
tau lepton at rest, Na is Avogadro's number and A is the 
nucleon number. The flux of muons depends on zenith 
angle 6 as well, but since we will look specifically at an- 
gles of = 0, 180°, we will not include the 9 dependence 
explicitly. A detailed parameterization of the 9 depen- 
dence appears in, for example, Ref. [llj. 

The differential energy distribution of the tau produc- 
tion cross section for Z^r is evaluated numerically follow- 
ing Ref. jJJl]. A parameterization for the decay distri- 
bution in terms of y = E^^ / E^ for Ztu^ can be found in 
Ref. 0. 

The atmospheric muon flux which appears in the /i — > 
T equation can also be represented by a transport equa- 
tion, but to first approximation, it can be written as 



<i>,,{E^,X)^c^^(ElQ)eMPX) . 
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In the above equation, 0^t(£'^,O) represents the muon 
flux at the surface of the Earth. The conventional at- 
mospheric muon flux from pion and kaon decay can be 
parameterized by jll| 
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Here 103 GeV and 810 GeV represent the pion and kaon 
critical energies, which separate the high and low energy 
contributions to the muon flux. The spectral index is 
7 = 2.72. 

We also include the muon flux from charm produc- 
tion and decay in the atmosphere, the so-called "prompt" 
flux. Here, we use the flux 
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from a dipole model evaluation of the charm production 
cross section [ij], although other choices are possible [J, 
0, 0]- The spectral index for the prompt flux is 7' = 
2.53. Note that the prompt muon neutrino flux equals 
the prompt muon flux. 

Taking energy losses through ionization and radiative 
processes into account, the muon energy at depth is re- 
lated to the surface energy through 

i?; = exp(/?X)£;^ + (exp(/3X)-l)| (8) 



in terms of the a — 2.67 x 10^ 
10-6 cmVg for E^ < 3.53 x 



'^ GeV cm^/g and /3 
10^ GeV, and a = 



-- 2 Ax 
-6.5 X 



lO-'^GeVcm^/g and /3 = 3.66 x 10~6cm2/g for higher 
energies [ll|- Putting Eqs. © and © back into Eq. © 
yields an expression for the muon flux as a function of 
energy and depth. 
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FIG. 1: Tau production Zfj,T moment for pair produced taus 
from incident muons. Shown are the Z,it moments for the 
conventional (solid) and prompt (dashed) muon fluxes with 
the approximations of Eqs. ^ and pop . 

An analytic solution to the cascade equations can be 
found when the fluxes can be factorized in the form 
(j)fj,{E,X) = (j)i{E)(j)2{X) 0. For the tau neutrino en- 
ergies of interest here, E^^ > 100 GeV, E^ > 100 
GeV, so the surface energy can be well approximated 
by £"* c± Ef^exp{(3X). For the depth ranges of impor- 
tance in our calculation, exp(— /JAT) ~ 1. Using these 
approximations, the incident muon fluxes can be factor- 
ized according to 
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The muon fluxes go into the evaluation of Z^r in Eq. 
([5]) . Figure [1] shows the tau production Z^^r moments 
using the approximate muon flux in Eq. ^ for the con- 
ventional (solid line) flux and Eq. (|TOl) for the prompt 
(dashed line) flux. The prompt Z-moment is larger than 
the conventional Z-moment because the prompt muon 
flux has a smaller spectral index. 



Using the factorized form of the muon flux at column 
depth X, the tau flux from the conventional muon flux 
as a function of energy and depth is 



-{Er,X) = 
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A similar expression is obtained for the tau flux from 
prompt muons. Because Xdec is a function of i?^, there 
is not a factorizable form applicable to all energies, 
(pr{E,X) ^ (f)i{E)(j)2{X) , so we have solved the cas- 
cade equation for the tau neutrino flux generated by tau 
decays numerically. 



III. OSCILLATED AND PROMPT 
ATMOSPHERIC TAU NEUTRINOS 

In addition to tau pair production and decay, there 
are two other sources of tau neutrinos. At low energies, 
the main source of tau neutrinos is from oscillations of 
i^fj. -^ Vr 14]. In the absence of matter effects, a good 
approximation for the energy range of interest [l5], the 
probability of i^^ — )■ ^'^- can be expressed in terms of the 
atmospheric mixing angle ^atm and mass splitting Ato^ 
in eV^ in the two flavor approximation. 



P{v^ -> Vr) = sin^(26'atm)sin^(l-27ATO 



E 



(12) 



For a distance L in km and energy E in GeV. For the 
atmospheric mixing angle and mass splitting, we use 



sin''(26'atm) = 1 and Am^ = 2.5 x lO"'^ eV 
calculate the tau neutrinos produced from oscillations of 
v^ — ^ Vt for the downward atmospheric flux, one needs 
to account for the production height. Gaisser and Stanev 
[31 have shown that for E^, > 100 GeV, the pathlength 
distribution (in £) of muon neutrinos is approximately 
flat out to icos9 ~ 20 km, followed by a declining pro- 
duction rate. We use the Bartol flux of muon neutrinos 
at sea level (j)^ = 4>i' +p {Ev,^) UI with a distribution 
that is constant up to L„iax cos 9 = 20 km. 
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In the high energy approximation with our choice of mix- 
ing angle and mass splitting, this gives a tau neutrino flux 
at the surface of the Earth of 
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where the factor of 1/3 accounts for the distribution in 
production heights. For a detector at depth d, looking 
up at a zenith angle 0, the oscillated tau neutrino flux is 
approximately 




X 10" 
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For our results of downward neutrinos, we consider a 
depth d ~ 1.5 km. 

For upward neutrino fluxes, the height of production 
and d are negligible compared to the length of the chord 
D through the Earth that the neutrinos travel, so we take 
the upward oscillated tau neutrino flux to be 






(16) 



The other source of atmospheric tau neutrinos is the 
decay of Ds mesons that are produced in cosmic ray inter- 
actions in the atmosphere [6|, |7| . An approximate form 
for the prompt flux of tau neutrinos and antineutrinos 
can be written [3] 

rr. r.^ Ix 10-7^°-5 (GeVcm^srs)"^ 
(t>u^+uAE,Q) = ,^;L.A. — (17) 
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At the energies of interest, the flux of prompt tau neu- 
trinos is independent of zenith angle due to the short 
lifetime of the I?, meson. 



IV. RESULTS AND DISCUSSION 

Our results for downward and upward tau neutrinos 
are shown in Figs. 2(a) and 2(b). Figure 2(a) shows 
the fi -^ T ^ Vr contributions of the conventional muon 
flux (solid line) and the prompt muon flux (dotted line). 



iM- lo in comparison to the oscillated conventional Vf^ 



flux for downward neutrinos with a zenith angle 9 = 0° 
(dashed line). With the short distance over which muon 
neutrinos can oscillate to tau neutrinos, muon production 
of taus which then decay, dominate over the oscillated 
conventional tau neutrino flux above E,y ~ 500 GeV. 

Figure 2(b) shows the upward {9 = 180°) tau neu- 
trino fluxes. There, with the oscillation distance equal 
to the diameter of the Earth, the oscillated flux domi- 
nates the /i — )■ r -H- J/,- over a much larger energy range. 
For upward neutrino fluxes, at the energies of interest, 

(13) neutrino attenuation needs to be considered. To account 
for neutrino attenuation of upward going neutrinos, we 
have included a factor of exp{—accXNA/A), where ace is 
the tau neutrino charged current cross section. For the 
charged current cross section, we have used the results 
from Ref. [l7| . This factor is a good approximation even 
for the /i ^^ r — > z^'t- flux because most of the tau pair pro- 
duction and subsequent decays happens near the surface 
of the earth. 

The dot-dashed lines in Figs. 2(a) and 2(b) shows the 
prompt i^T + S't flux evaluated in Ref. [7|] given by Eq. 
(fTT)) from Ds production followed by its decay to rvr- 
While there is some uncertainty in the prompt tau neu- 

n 5") trino flux, the figures show that it dominates the atmo- 
spheric contributions to the flux of Vr + v^- over a wide 
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range of energies discussed here. The dominant contribu- 
tion to the tau neutrino flux for i? > 2 TeV [1^1, regard- 
less of the zenith angle direction, is the direct prompt pro- 
duction of tau neutrinos. For downward neutrinos, the 
direct prompt production dominates the other sources of 
tau neutrinos for energies larger than tens of GeV. 

In addition to downward {9 = 0°) and upward {0 = 
180°) tau neutrino fluxes, the flux of tau neutrinos at 
higher zenith angles can be calculated. For higher zenith 
angles, the conventional fluxes are increased, while the 
prompt fluxes are essentially unchanged for .E < 10^ 
GeV. The increased conventional muon flux means that 
1000 10* 10° lo'' there is an increase in the flux of neutrinos produced via 

El, [GeV] ^ ^ Y ^>- Ur- Similarly, however, there is an increase 

in the flux of tau neutrinos from the longer pathlength 
for oscillations plus the higher flux of muon neutrinos. 
At an incident zenith angle of 80°, tau neutrinos from 
/i — >■ r — >■ i^T- dominate those from oscillations over the 
approximately the same energy range as for 6 — 0°. At 
9 = 80°, the conventional fi ^- t -^ Vr and i'^ -^ i>r 
oscillated fluxes are increased by a factor of about three, 
but the prompt i^r flux still dominates. 

In summary, while the flux of downward tau neutrinos 
produced from atmospheric muons dominates that pro- 
duced by neutrino oscillations for a large range of energies 
and zenith angles, both of these sources are suppressed in 
1000 10* 10° 10" comparison to the flux of prompt tau neutrinos coming 

E^ [GeV] from the decay of Dg mesons. 



(b) 



FIG. 2: Differential tau neutrino flux scaled by the square of 
the neutrino energy for neutrinos (a) going downward (with a 
zenith angle of 0°) and (b) going upward (with a zenith angle 
of 180°). The solid curve represents the tau neutrino flux 
from the incident conventional muon flux (/i — )■ r — s- !/t), and 
the dotted curve comes from the prompt atmospheric muons. 
The dashed curve shows the oscillated conventional flux from 
muon neutrinos and the dot-dashed curve shows the prompt 
tau neutrino flux itself. 
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